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ABSTRACT 

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is 
frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of 
this approach analytically and numerically. We argue that the character of the flow 
where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle 
model. Only in the specific case of hot virialised gas with zero angular momentum 
and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite 
extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient 
cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism 
can be erroneous by orders of magnitude in either direction. This may impose artificial 
trends with halo mass in cosmological simulations by being wrong by different factors 
for different halo masses. We propose an expression for the sub-grid accretion rate 
which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore 
taking account of the contribution of the halo to the gas dynamics. 
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1 INTRODUCTION 

Over the last decade, compelling observational evidence 
has revealed that many galaxies in the local Universe 
harbour supermassive black holes (SMBHs) with masses 
10^ Mbh/M© ^ 10^ in their centres. During the same pe- 
riod, surveys of the distant Universe uncovered the existence 
of quasars at ^ ~ 7, when the Universe ^ I/IO"^ of its cur- 
rent age; this imphes that ma ny SMBHs had alread y assem- 
bled their mass by this time (jMortlock et al.. 201ll ). 

Our understanding of the physics that dictates the 
growth of SMBHs is incomplete. Black holes grow by accret- 
ing low angular momentum material from their surround- 
ings, yet the character of the accretion flow onto an SMBH 
is gover ned by physical processes as diverse as galaxy merg- 
ers (e.g.. lHopkins fc Quataert. 20101'). turb ulence induced by 
stellar feedback (e.g.. iHobbs et al.. 201l|) and black ho le 
accretion-driven outflows (e.g., iNavakshin fc Power. 2O10l '). 

Black hole growth is now routinely modelled in 
galaxy formation s imulations (for a flducial work see 
ISpringel et al.. 20051 ) and the importance of SMBHs 
in shaping the pr operties of galaxies is now well es - 
tablished (e.g. ICroton et al.. 200dBower et al.. 20061 ). 
The majority of galaxy formation simulations pub- 
lished in the literature incorporate what we shall term 
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the "Bondi-Hoyle model" for black hole growth (see, e.g., 

ISpringel et al. . 200dSiiacki et al. 200iPelupessy et al 2007^ Matteo et al. 



which derive s from the work of iBondi 
iBondi (1951 ) 



Hwle (19441 ) and 
This model assumes the 



hereafter B&H. I'liis model assumes 
simplest possible accretion flow, where the gas is at rest at 
inflnity and accretes steadily onto a black hole, subject only 
to the (Newtonian) gravity of the latter, which is modelled 
as a point mass. 

Simulations that model the idealised physical problem 
as it is set out in B&H, or in idealised generalisations, 
produce results that are in good agreement with the 
analytical solution (|Ruffert. 1994Barai et al.. 201ll ). In 
galaxy formation simulations, unfortunately, this idealised 
picture is far from satisfled, as the gas inflow is compli- 
cated considerably by the properties of the flow at larger 
scales. The most notable example of this is the presence 
of non-zero angular momentum that provides a natural 
barrier to eventual accretion by the SMBH. Gas settles into 
a disc whose dimensions are set by the angular momentum 
of the accretion flow, with only the very lowest angular 
momentum material able to accrete. A true estimate of 
the accretion rate onto the SMBH must therefore take 
account of this angular momentum, and indeed attempts to 
include it in an accretion sub-grid model have been ma de 
iPower et al.. 20ldLevine et al. 201CtDeBuhr et al.. 201 j ) 
along with t he presence of turbulence and /or vorticit; 
in the gas jKrumholz et al.. 200dKrumholz et al.. 200 
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finding large departures from the standard Bondi-Hoyle 
rate. 

In this short paper we wish to make a simple and more 
fundamental point that in galax;y formation simulations even 
spherically symmetric accretion cannot be correctly mod- 
elled by the Bondi-Hoyle formalism, except in the most spe- 
cific of cases. To make this point we suspend, for the mo- 
ment, our disbelief that gaseous infall can proceed entirely 
radially from large scales and consider zero angular momen- 
tum accretion fiows onto an SMBH embedded in the poten- 
tial of a massive dark matter halo. Indeed, this is an example 
of a situation where one might expect the B&H formula to 
provide a reasonable estimate of the accretion rate. 

The layout of this paper is as follows. In Section 2 we 
show analytically that in large-scale simulations of cosmo- 
logical volumes the Bondi-Hoyle approach is invalid, and in 
Section 3 we present some numerical tests of this hypothesis. 
Finally in Sections 4 & 5 we discuss our conclusions. 



2 ANALYTICAL ARGUMENTS 

2.1 Classical Bondi-Hoyle accretion 

We first recap the main assumptions underpinning the clas- 
sical B&H pape rs. These ar e nicely summarised in the first 
sentence of the lBondi fl952l Vs abstract: "The special accre- 
tion problem is investigated in which the motion is steady 
and spherically symmetrical, the gas bemg at rest at infin- 
ity" . We have italicised the part of the sentence that bears 
the most importance for us here. 

Physically, gas can be at rest at infinity only when it is 
not subject to any forces. The only external force acting on 
the gas in the restricted B&H problem, i.e., the gravitational 
force, is due to the black hole. Self-gravity of the gas is 
neglected. The "infinity" in question is a region at a distance 
large enough from the SMBH that the gravitational force 
exerted by the latter is negligible when compared to the 
pressure forces within the gas. This is quantified by defining 
the Bondi (or the "capture") radius. 



Applying the Bondi-Hoyle formalism to black hole 
growth assumes that the accretion rate onto the SMBH is 
commensurate with the accretion rate through the Bondi 
radius (i.e., the flow is steady-state) and therefore given by 



TB = 2GMbh/4 



(1) 



where Mbh is the mass of the central object and Coo 
is the sound speed of the gas far from the hole. The 
Bondi radius divi des the flow into two distinct regions 
jFrank et al.. 20021 ). Far from rs, gas is hardly aware of the 
existence of the black hole, and the flow is very subsonic. The 
pressure and density of a subsonic flow are approximately 
constant, therefore we can set p(r) « paa at r ^ rs. 

Inside the capture radius, on the other hand, p(r) be- 
gins to increase above the initial value, and the flow even- 
tually reaches a sonic point where \vr\ — Coo, within which 
it plunges essentially at free-fall. The sonic point is found 
from Ts = GMBH/2Cs(rs), where Cs(rs) is the sound speed 
at Ta. T his local quantit>n is related to the sound speed at 
infinity l|Frank et al.. 2002 ) via Cs{rs) = Coo(2/(5 - ar))^''^ 
where F is the polytropic index of the gas that relates the 
gas pressure and density by P = Kp^ , with K a positive 
constant. 



note that r's here is the sonic radius, not the scale radius as it 
is commonly used in descriptions of dark matter halo profiles 



Mbh = 7rA(r)rBPooCo 



47rA(r)G2MiHpo 



(2) 



where A(r) contains all the corrections arising due to the 
finite pressure gradient force in the problem. This function 
varies relativ ely weakly, i.e ., between 1.12 for F = 1 and 0.25 
for F = 5/3 iBondi. 19521 ). For the remainder of this paper 
our fiducial assumption is a soft equation of state i.e., F « 1. 



2.2 When is Bondi-Hoyle accretion applicable? 

The Bondi-Hoyle formalism has been widely adopted as a 
'sub-grid' prescription for the accretion rate onto the SMBH 
in large-scale cosmological simulations. The standard argu- 
ment is that while one cannot usually resolve the scales 
of the Bondi radius, one can at least use the smallest re- 
solved scales to approximately determine the value of the 
gas density and the sou nd speed "at infinity" f or use in the 
Bondi formula (see, e.g.. lBooth fc Schave. 20091 ). Indeed, the 
smallest resolvable scales are usually about a fraction of a 
kpc, whereas the Bondi radius is of the order of a few to a 

few tens of pc^ 

For the lBondi (19521 ) solution to be applicable even to 
spherical flow, we need to make sure that the "gas being 
at rest at infinity" assumption is satisfied where the rel- 
evant gas properties (density and sound speed) are eval- 
uated. In cosmological simulations we expect SMBHs to 
be immersed in stellar bulges and dark matter haloes that 
are typically ~ 1 0'^ to lO'* times more massive than the 
SMBH (see, e.g., iHaring fc Rix. 2004Guo et al.. 2010| ). If 
gas in the halo (or bulge) is as hot as the halo virial 
temperature, it will be in hydrostatic balance (see, e.g., 
iKomatsu & Seliak. 200llSuto et al.. 1998t ). It seems that 
in this situation, which is common for low luminos- 
ity SMBHs in giant elliptical galaxies where the gas is 
rather tenuous and h ot since the cooling time is long 
(|Churazov et al.. 2005I ) the Bondi-Hoyle solution is poten- 
tially useful. 

However, in the epoch when SMBHs grow rapidly, 
their hosts are very gas rich, and the inflow of 
gas from large scales cannot be easily captured by 

gradual cooling from a tenuous hot halo (see, e.g., 

iBirnboim fc Dekel. 200dKeres et al.. 200dKeres et al.. 200gDekel et al.. 200^ 

Higher density gas is likely to cool much faster and hence 

is likely to be much cooler than the virial temperature. In 

this case the gas is not able to support its own weight, and 

must collapse to the centre, where it feeds the SMBH and 

forms stars. Therefore we expect a radial inflow of gas to 

the centre rather than a hydrostatic b alance "far" from the 

SMBH. We note that iRicotti (20071 ) have demonstrated 

how the Bondi-Hoyle formalism must be modified in the 

presence of an external potential for the specific case of the 

growth of primordial black holes in a dark matter halo with 

a power-law density profile. 

In actuality, even the meaning of the Bondi radius be- 
comes unclear in this situation, as gas is not virialised near 
the SMBH, and the halo potential plays an important role. 
In the "naked SMBH" problem, rs delineates the region 
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inside which the potential energy of the hole starts to be- 
come greater than the internal energy of the gas. For an 
SMBH plus host halo system, one should introduce a mod- 
ified Bondi radius, 

fs = 2GMn/ci, (3) 

that takes into account the total mass of the halo. In or- 
der for the gas to accrete efficiently onto a dark mat- 
ter halo, its temperature must be at most comparable 
with the virial temper ature at the outer edge of the halo 
(|White fc Frenk. 1993 ). Thus, we set Coo ~ GMh/vh, where 
Th is the halo virial radius, giving us a modified Bondi ra- 
dius of fs XL 2rh- The gravitational potential energy starts 
to dominate the internal energy of the infalling gas before 
the latter has even reached the edge of the halo, and so the 
standard rs is meaningless in this case. 

Figure [T] illustrates this point graphically by comparing 
the potential energy of g function of radius for a 

variety of halo profiles with assuming that the gas 

temperature is virial at r^. These profiles are described in 
the Appendix. The cosmology we have assumed for the halos 
is ACDM with flm = 0.27, fl\ — 0.73 and a virial overdensity 
parameter of A = 200 at a redshift of z — 2 (although we 
note that our conclusions are unchanged for a wide range in 
z, from the early Universe to the present day). 

The oc 1/r potential due to the SMBH of mass Mbh = 
10* M0 is shown with the dash-triple-dot power law. The 
traditional Bondi (or capture) radius is at r = 0.002, where 
the SMBH potential crosses the c5ir/2 line. For all the halos 
considered, the potential energy of the host makes a signif- 
icant contribution at all radii except those within the very 
inner parts of the halo, exceeding the gas internal energy 
by a large factor everywhere outside the SMBH radius of 
infiuence. 

Therefore, if we assume that the gas has a relatively 
soft equation of state and accretes spherically onto the halo 
at or below Tviiiai{rh), the character of the infiow becomes 
that of supersonic inside the halo instead of bemg stationary 
at infinity. In the Appendix we quantify this by calculating 
the sonic point for isothermal gas flows at or below the virial 
temperature and for all of the halo mass profiles considered. 
We show that for typical dark matter halos the sonic point 
is reached while still at very large distances from the central 
black hole. 

From Figure[T]it can also be seen that for hotter gas, the 
classical tb estimate becomes more accurate, as the poten- 
tial energy of the (halo + SMBH) system starts to asymp- 
tote to the SMBH solution and no longer dominates over 
the value of cJii./2. If the gas thermal energy is comparable 
to that of the potential energy we naturally expect near hy- 
drostatic equilibrium to be maintained and in this case the 
Bondi-Hoyle formalism is applicable to spherical fiow. 

Of course, this is a simplified picture. In reality gas 
that has accreted onto a halo from outside the virial ra- 
dius may shock at smaller radii, heating up to the local Tvir 
at that radius. In the 'cold mode' where the galaxy is as- 
sembled via cold streams that penetrate far inside the halo, 
such shock heating is likely to occur when the infalling gas 
reache s the radius of the galactic disc, at a small fraction 
of iKeres et al.. 20051 ). However, due to the high den- 
sities reach ed for the shocked gas the cooling time is likely 
to be short (|Birnboim fc Dekel. 200dKeres et al.. 20051 ') , and 




rB@T=1.4x10°K 



1.e-05 l.e-04 1.e-03 1.e-02 1.e-01 1.e+00 1.e+01 

Figure 1. Plot of potential energies for a selection of dark mat- 
ter halo profiles (some realistic, some instructive) modified by a 
central SMBH, compared to the thermal energy of the gas at the 
virial temperature at the halo virial radius. The halo profiles are 
normalised to contain a mass = 10^^ Mq within the halo 
radius r^i, outside of which the potential is Keplerian by con- 
struction. The potential due to the central SMBH by itself is also 
plotted. The classical Bondi radius for three different values of 
temperature are indicated as crosses - the greater the separation 
between the cross denoting rs and the curves of potential energy 
for a given halo the more inaccurate the standard Bondi-Hoyle 
approach. It is clear then that the latter becomes more accurate 
at higher T. The modified Bondi radius for the fiducial T^i^ is 
denoted by a rhombus, demonstrating that at this temperature 
the potential energies are dominant for the entirety of the halo. 



in particular in the presence of Compton cooling from a 
quasar radiation field it will be significantly less tha n the 
free-fall time at that radius (jNulsen fc Fabian. 200G| ). Gas 
may therefore be stationary temporarily at the shock radius 
but as it cools and begins to infall it will quickly tend to the 
free-fall velocity, and certainly by the time a radial infiow 
reaches the classical Bondi radius - values for which can be 
seen in Figure [1] for different temperatures - the assumption 
of being at rest will no longer be satisfied. 

Put simply, this is an energy argument. Regardless of 
where the gas might begin to infall from, the halo imposes a 
far stronger gravitational potential energy than the SMBH 
everywhere outside of the SMBH radius of influence (by def- 
inition). The potential energy due to a typical halo also in- 
creases to smaller radii from any point at which the gas is 
likely to have reached the virial temperature (see Figure [l|. 
As a result it is clear that gravity will often dominate over 
thermal energy at scales significantly larger than the classi- 
cal Bondi radius for the SMBH. 

In the interest of completeness we now demonstrate nu- 
merically the form of the radial infall in a realistic back- 
ground potential within a galaxy. We choose a dynamic 
range that lies inside the sonic point (as would always be 
the case for efficiently cooled gas at these scales, as we have 
shown in the Appendix) in order to highlight why the stan- 
dard Bondi estimate is inaccurate here. 
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3 NUMERICAL TESTS 

To perform the simulations we employ the three-dimensional 
smoothed particle hydrodynamics (SPH)/N-body code 
GADGET-3, a n updated version of the code presented in 
ISpringel f20"05l ). The gas is evolved in a static external po- 
tential that includes a point mass black hole at the centre. 
The computational domain extends from a kiloparsec down 
to an "accretion radius" around the black hole at race ~ 1 
pc, and we remove the particles that come within this dis- 
tance of the SMBH. 

For the external potential in our model we use a Jaffe 
cusp as per equation I A3I but with a core at the centre of our 
computational domain to prevent divergence in the gravita- 
tional force. The radius of the core, rc, corresponds to ap- 
proximately the dynamical influence radius of the SMBH. 
With this (modified) potential the mass enclosed within ra- 
dius r is given by: 

M(r) = Mm.H- { ' . (4) 

where Mc = 2 x 10** Mq, M = lO" Mq, Tc = 20 pc, and 
a = 10 kpc. The mass of the SMBH is set to Mbh = 10** Mq. 

For simplicity, the gas is kept isothermal throughout the 
entirety of the simulation and self-gravity is turned off. 





' 1 1 ' ' ' 

(B-H)T=1e3 
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Eq. 5 
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Figure 2. The mass flux (steady-state value) from the simu- 
lations through each radius (red) i.e., M(r), together with the 
Bondi-Hoyle estimate as per Equation[2]and the free-fall estimate 
as per Equations [5] & [6] 

parts results in a significant overestimate of the accretion 
rate and an underestimate at large radii. 



3.1 Initial conditions 

The starting condition for our simulations is that of a uni- 
form density, spherically symmetric thick gaseous shell with 
mass Mahcii = 10® Mq, that ranges from rin = 0.1 kpc to 
Tout = 1 kpc and is centered on the black hole. The temper- 
ature of the gas is varied between tests, ranging from 10^ 
K to 10^ K. To minimise initial inhomogeneities we cut the 
shell from a relaxed, glass-like configuration. 

The gas infalls from rest within our static potential. 
After a time of the order of the dynamical time at the outer 
edge of the shell, a steady-state radial mass flux is reached 
for the majority of the gas and it is at this time that we 
make our comparisons between the various accretion rates 
in the next section. 



3.2 Results 

We define a radius-dependent "measured" accretion rate as 
M(r) = 4-Kr^pVr, where Vr is the radial velocity of the gas. 
Since the system has been allowed to settle into an approxi- 
mate steady-state, this function (shown as an average by the 
long dashed red line) is almost constant with radius, and is 
the same as the time-averaged accretion rate measured at 
the black hole. 

The solid, the dotted and the dashed curves in Figure 
[2] show the standard Bondi-Hoyle estimate for the accretion 
rate (equation [2| as a function of radius for three different 
values of gas temperature. Each of these curves uses a 'local' 
p, the density at each radius, in order to represent where the 
poo might be evaluated in a large-scale simulation. 

Clearly, as Figure [5] shows, the Bondi-Hoyle estimate is 
very inaccurate for these isothermal simulations. At inter- 
mediate temperatures, e.g., 10^ K, the formula in the inner 



3.2.1 Free-fall rate 

A simple but physically well motivated alternative to 
the Bondi-Hoyle formula for efficiently cooled spherically- 
symmetric flows is a free-fall rate estimate, 

. ^l^SgpM , (5) 

where Mgas,enc(r) is the enclosed gas mass within radius r, 
and tft = {r^ /2GM{r)y^^ is the free faU time. Alternatively 
one may wish to approximate the enclosed gas mass in the 
above equation as (47r/3)r^/9gas(r'), so that 

Both of the above estimates provide a much better 
match to the accretion rate in our numerical tests than the 
Bondi-Hoyle approach, as Figure [2] shows. It should be noted 
too that since the free-fall estimates have no dependance on 
Ca, the profiles are converged regardless of the temperature. 



4 DISCUSSION 

Based on our anal ytical argume nts and numerical tests, we 
conclude that the iBondi fl952l ) formula, designed for ac- 
cretion onto "naked black holes", can only be applied to 
accretion onto astrophysical SMBHs, i.e., those embedded 
in massive dark matter halos, if the gas in the latter is at 
or near hydrostatic equilibriu m. It is only i n this case that 
one of the key assumptions of lBondi (19521 ) - the "gas be- 
ing at rest at infinity" - is satisfied (infinity meaning outside 
the Bondi radius). If the gas cooling time is long then it is 
possible that this state may be reached, and indeed this is 
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probably the case in giant gas-poor elliptical galaxies, wher e 
the gas is tenuous and hot (see, e.g. JChurazov et al.. 2005h . 

However, in the most interesting phase of SMBH and 
galaxy buildup, when the halo is likely to be awash with gas 
to feed both the SMBH and star formation, densities are 
high and the cooling time is expected to be short, with gas 
cooling far below the virial temperature. Hydrostatic bal- 
ance is then extremely unlikely for gas in the halo, and if it 
is not supported by either shock heating or angular momen- 
tum it will tend to free-fall velocities. We have run a series 
of simple numerical tests to explore this limit, allowing a 
thick, spherical shell of gas to accrete onto an SMBH at the 
centre of a background halo potential. These tests showed 
that a free-fall accretion rate estimate is indeed much more 
accurate than the Bondi-Hoyle formalism in this case. What 
is most concerning is that the error of the latter strongly de- 
pends on gas temperature (and thus the cooling function) 
and cannot be "predicted". The Bondi-Hoyle accretion for- 
malism may thus be wrong by significant (and unknown) 
factors in either direction. The error may well depend on 
the galaxy mass and type systematically. Predictions based 
on the Bondi-Hoyle formalism alone are therefore unlikely 
to be robust. 

We note too that this problem persists even for simu- 
lations that are able to resolve down to the classical Bondi 
radius of the SMBH. While the properties of the gas (den- 
sity, sound speed) evaluated at this radius may be correct in 
the sense that they do not suffer from resolution problems, 
it is unlikely that this radius will be a true "infinity" where 
the gas is in hydrostatic equilibrium - again, the reason for 
this is the influence of the more massive background halo 
potential on the dynamics of the gas. 

For realistic flows of course the situation is more com- 
plex. Here we find that the best way to present our argument 
is in terms of "supporting" mechanisms for the gas at the ra- 
dius where one evaluates the gas properties in a simulation. 
To start with, we have the two extremes: (i) Bondi-Hoyle, 
where the gas is completely supported and stationary at the 
evaluation radius; and (ii) free-fall, where the gas is entirely 
radial and influenced only by gravity i.e., has zero support. It 
is important to recognise that the reality will lie somewhere 
in-between. Specific cases that lie between these extremes 
are: (iii) the gas is unsupported at the evaluation radius but 
is shocked at smaller radii, after which it may begin to tend 
to free-fall once again if the cooling time is short, and (iv) 
the gas possesses angular momentum, and is thus supported 
by rotation but inflows through viscous processes. 

What we therefore desire is a formula for the sub- 
grid accretion rate that interpolates between the two 
extremes based on the relevant supporting mechanism. 
As we mentioned in the introduction, angular momentum 
concerns are extremely important in determining the 
correct accretion rate; and indeed the picture can be further 
complicated by the effect of star formation, where it is not 
clear whether forming stars would deprive the SMBH of 
fuel (as may have been th e case with Sgr A*- see, e.g., 
iNavakshin fc Cuadra. 20051 ) . or if feedback from star forma- 
tion would actually enhance accretion t hrough a broadenin g 
of the angular momentum distribution (|Hobbs et al.. 2011f ). 
Feedback from the AGN itself must also be a consideration, 
and in particular the interplay between SMBH feeding 
and AGN feedback, due to the fact that these processes 



connect the small to the large scales. Indeed, a number 
of authors have conducted detailed investigations into the 
multi-scale coupling between feeding and feedback (see, e.g., 
ICattaneo fc Tevssier. 200iDubois et al.. 201(lKim et al.. 20111 ), 
often finding that SMBH growth enters into a self-regulated 
state that oscillates between periods of high and low activ- 
ity. We plan to investigate all of the processes relevant to 
SMBH growth and thereby to further develop the sub-grid 
model; for now, however, we focus on a direct interpolation 
between the Bondi-Hoyle regime and the free-fall regime as 
a starting point. 

The formula we propose is a minor mod- 
ification of the full Bondi-Hoyle-Lyttleton 

(|Hovle fc Lvttleton. 193dBondi. 1952h expression for 
an accretor that is moving relative to the gas, namely 

47rA(r)G'"MiHpoo 



Mbhl = 



(7) 



where v^ai is the relative velocity between the SMBH and 
the gas, which in purely radial flow is zero. In order to 
interpolate between this and the free-fall rate (equation 
O we make two changes: (1) replace the relative veloc- 
ity with the velocity dispersion for the external potential, 
(T ~ (GMenc{r) /ry^'^ , and (2) replace the black hole mass 
with the enclosed mass of the external potential, Menc(^')- 
The resulting expression. 



(cL+a2)3/2 



(8) 



tends to the standard Bondi-Hoyle formula (equation [21) in 
the limit that c^o 3> cr, i.e., the gas internal energy dom- 
inates over the potential energy of the halo, and Mcnc — >■ 
Mbh, i.e., as we approach the SMBH radius of influence. In 
the opposite limit, where the halo potential energy domi- 
nates, a 3> Coo, we recover the free-fall estimate, equation 
[S] Finally, in order to incorporate (in a crude fashion) the 
effect of AGN feedback on the accretion the rate should be 
capped at Eddington, namely 

Afacc = max(Mintorp, MEdd) (9) 

as is commonly employed in the majority of the simulations 
of SMBH growth that we have mentioned so far. 

We hope that with this interpolated expression the sub- 
grid SMBH accretion rate will take account of the presence 
of the halo and indeed any other external potential (e.g. stel- 
lar bulge, stellar halo, etc.) that contributes signiflcantly to 
the dynamics of the infalling gas. We note that some simu- 
lations may already use a velocity dispersion as a cap on the 
v^ci parameter in order to avoid exc essive relative veloc ities 
as a result of poor resolution (e.g.. [Dubois et al.. 2010| ). in 
which case the change to the sub-grid expression is minimal. 

This interpolation approach is particularly relevant 
to the picture of the 'hot' and 'cold' modes of accretion 
in galaxy formation, that emerged in the 1970s with the 
paradigm of hot halo gas that is at or near Tvir after being 
shock heated at the virial radius (see the original papers b y 
e.g., iRees fc Ostriker. 197iSilk. 197lWhite fc Rees. 19781 ) . 
and has since been adjusted by results from simulations 
(although was suggested also in the 19 70s based on an- 
alytical arguments - see iBinnev. IQTtI ) to include cold 
gas with a soft equation of state that penetrates down 

to small scales through streams and filaments (see, e.g., 

iKatz et al.. 1994Fardal et al.. 200]|Keres et al.. 200dKeres et al.. 20091 ). 



© 0000 RAS, MNRAS 000, 000-000 



6 Alexander Hobbs, Chris Power, Sergei Nayakshin & Andrew R. King 



There is typically a halo mass scale which delineates the 
relative importance of each mode, although this varies some- 
what in the literature, with iBirnboim fc Dekel (20031 ) find- 
ing Mhaio '--^ lO'^^ as the result of 1-D numerical calculations 
and iKeres et al. (2005h finding a factor of 2-3 higher from 
fully 3-D N-body/hy drodynamical simu lations. However, in 
a subsequent paper, iKeres et al. iiooi) find that although 
this transition mass marks the point above which hot, viri- 
alised gas atmospheres develop in halos, the actual contri- 
bution to the accretion rate to small scales is still dominated 
by the cold filamentary mode even at higher mass, with hot 
mode accretion only starting to become important at late 
times (z < 1). In addition, there is (indirect) observational 
evidence for cold accretion flo ws in galaxies through th e 
HI column density distribution (jvan de Voort et al.. 201lh . 
This suggests that utilising purely a 'gas being fully sup- 
ported' Bondi-Hoyle accretion prescription is unlikely to 
capture enough of the relevant accretion behaviour in galaxy 
formation. Using the interpolation expression above would 
automatically adjust for the hot and cold mode dominance, 
via the relative importance of the and terms in 
the denominator. We note that semi-analytic models of 
galaxy formation often employ a distinction in the treat- 
ment of accretion rate between the two modes (see, e.g., 
iHirschmann et al.. 201lh . 

We would like, however, to emphasise again the most 
important caveat with the picture we have presented - the 
lack of angular momentum. If indeed the free-fall mode of 
accretion becomes more dominant as a result of employing 
the interpolated accretion rate expression then one must 
be very careful, for it is possible that the residual angu- 
lar momentum of gas infalling at free-fall velocities from 
large scales becomes even m ore relevant to the accretion 
rate (|Nulsen fc Fabian. 200(il '). due to the formation of a cen- 
trifugally supported disc at small scales. In this case we 
argue that any expression which does not take account of 
angular momentum should be viewed as more of a 'cap- 
ture' rate, with the actual accretion onto the SMBH mod- 
elled with a viscous timescale, the parameters of which are 
set by the properties of the gas at the transition point be- 
tween 'infai r and 'disc-mode' provided it is sufficiently re- 
solved (see iPower et al.. 2010| '|. We note for completeness 
however that an alternative picture where the formation of 
a disc is not in fact a hinderan ce to SMBH growth was 
presented by iMaver et al. (20101 ') ■ who found that global 
instabilities in the disc-like structure that formed from a 
merger led it to collapse and feed the SMBH on a dynam- 
ical timescale (and see als o the multi-scale simulations by 
iHopkins fc Quataert. 2010l ). 



5 CONCLUSION 

In this paper we have shown that the de-facto industry stan- 
dard - the Bondi-Hoyle formalism for accretion rate onto the 
SMBH fails for more than one reason in a realistic cosmo- 
logical simulation whenever gas cooling is efficient. We have 
shown that a free-fall estimate is more appropriate in this 
case (provided that angular momentum does not impede ac- 
cretion of gas even further). We suggested an approximate 
interpolation formula that bridges the rapidly cooling and 
the inefficient cooling regimes which we hope will be use- 



ful for cosmological simulations that cannot resolve gas flow 
all the way down to the Bondi radius (a few to few tens of 
parsecs). 

Finally, we conclude with the answer to the question 
posed in the paper title, which we feel is worth re-iterating: 
what matters in a simulation is the character of the flow 
where one evaluates the accretion/capture rate. If it is (a) 
spherical, (b) at rest, and (c) sufficiently far away to count 
as 'infinity', with (d) the enclosed mass dominated by the 
SMBH, and finally (e) the flow is uninterrupted between the 
evaluation radius and the black hole, then and only then is 
Bondi-Hoyle accretion applicable as a sub-grid model. As we 
have mentioned, these latter requirements may be met to a 
sufficient degree by massive gas-poor ellipticals if sufficiently 
resolved in a simulation, but are not by any situation where 
there is appreciable infall from large scales. We note that 
the interpolation formula we have proposed is applicable to 
both cases, although requires further development in order 
to take account of the full range of accretion regimes as 
discussed. 



6 ACKNOWLEDGMENTS 

We thank the referee for a detailed report that helped im- 
prove the paper. We acknowledge useful discussions with 
Justin Read, and thank Volker Springel for the version of 
GADGET-3 that was used for the simulations. Theoretical 
Astrophysics at the University of Leicester is supported by 
an STFC rolling grant. 



APPENDIX A: THE SONIC POINT 

For an isothermal gas flow, the sonic point for an extended 
mass distribution satisfies r = GM{r)/2cl{r), where M(r-) 
is the enclosed mass. Depending on the mass profile, this 
may be inside or outside the halo. We now consider several 
profiles for dark matter halos, solve for the sonic point, and 
plot its location versus gas temperature in Figure lAll In 
particular, for a power-law density profile with index q, 



M{r) = Mh - 



■A-q 



(Al) 



SO that Ts = 2^^^"^ 'Vh for gas at rvir(j"h). The enclosed 
mass in the Navarr o, Frenk and White (NFW) profile 
(|Navarro et al.. 19961) is 



M(r) = 47r/5oa 



In 1 + 



a{r + a) 



(A2) 



where a is the scale radius of the halo and po is a characteris- 
tic density. The scale radius a depends on the concentration 
c = rh/a of the halo , where 5 < c < 15 for halos of mass 
10" -10^° Mo (e.g.. iBullock et al.. 200ll ). 

Elliptical galaxies an d bulges that follow the B}^^ law 
of Hi Vaucouleurs (1948t ) can usually be modelled by one 
of a family of spherical den sity profiles characterised by an 
exponent 7 (jDehnen. 1993h . Here the enclosed mass goes as 



M(r) = A-Kpoa 



r + a 



3-7 



(A3) 



with the Hernquist profile (|Hernquist. 1990l ) and the Jaffe 
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q = 3 power law 

q = 4 power law 

Hernquist 

NFW 

Jaffe 



1.0 

T/T„, (r,) 

Figure Al. Sonic radius, scaled in units of rj,, for a variety of 
halo profiles, assuming a halo of lO^'^ Mq. 

cusp corresponding to 7 = 1 and -y = 2 respec- 

tively. 

Referring now to Figure [Al] we can see that the sonic 
point is in all these cases a significant fraction of the halo 
radius at T = Tvir{rh), and becomes larger than rh for cooler 
gas. Thus, for a spherically-symmetric problem, the infalling 
gas at or below the virial temperature (defined at rh) must 
quickly tend to a free-fall solution once inside the halo. 
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